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Abstract — The work reported here details basic validation of 
our advanced physics-based EMI forward and inverse models 
against data collected by the NRL TEMTADS system. The data 
was collected under laboratory-type conditions using both artifi- 
cial spheroidal targets and real UXO. The artificial target models 
are essentially exact, and enable detailed comparison of theory 
and data in support of measurement platform characterization 
and target identification. Real UXO targets cannot be treated 
exactly, but it is demonstrated that quantitative comparisons of 
the data with the spheroid models nevertheless aids in extracting 
key target discrimination information, such as target geometry 
and hollow target shell thickness. 

I. Introduction 

CLEANUP of buried unexploded ordnance (UXO) from 
old practice ranges is a longstanding economic and 
humanitarian problem. Solution of the problem requires re- 
mote identification of subsurface metallic objects. The most 
difficult technological issue is not the detection of such 
targets — advanced metal detection systems, such as the NRL- 
TEMTADS platform described below, easily detect even very 
small amounts of metal at 1 m or more depths — but rather the 
ability to distinguish between targets of interest and harmless 
clutter items, such as various sized fragments of exploded 
ordnance. Since clutter tends to exist at much higher density, 
even modest discrimination ability leads to huge reductions in 
the economic cost of remediating such sites [1]. 

Formally, a successful solution to the electromagnetic (EM) 
discrimination problem is an algorithm enabling accurate 
bounds on physical properties of the target scatterer (position, 
shape, orientation, composition, etc.) from measurements of 
the scattered field using a well characterized platform (known 
transmitter/receiver geometry, transmitted waveform, and so 
on). Solution of this inverse problem requires a search over 
candidate solutions to the forward problem, namely accurate 
forms for the scattered field from a known target in a known 
subsurface environment. Generating high-fidelity forward so- 
lutions requires full three dimensional numerical solutions 
to the Maxwell equations, a difficult and time consuming 
computational problem. To reduce the computational burden, it 
is extremely important to obtain analytic solutions to as broad 
an array of exactly soluble model problems as possible. These 
solutions may then either be used as first-order models of the 
target, or as the basis of a perturbation scheme for accurate 
modeling of "nearby" target geometries. 

This paper details successful validation of our physics-based 
"mean field" and "early time" approaches to modeling of 
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time-domain electromagnetic (TDEM) responses of compact, 
highly conducting targets. Specifically, we apply our methods 
to the analysis (through both forward and inverse modeling) of 
laboratory-style data collected by the NRL TEMTADS system 
using artificial spheroidal targets, as well as some real UXO 
targets. The models use the detailed measurement platform and 
target parameters to generate highly numerically efficient, first 
principles predictions for the measured time-domain voltages. 
The models are designed to be essentially exact for spheroidal 
targets, and the remarkable agreement between measurements 
and predictions strongly supports this conclusion. The EM 
response of real UXO targets is found to differ in significant 
ways from those of spheroids, but comparing the two provides 
key insights into target identification. 

The only compact targets for which a full analytic solution 
at any frequency may be derived are those with spherical 
symmetry [2]. These are rather poor approximations to UXO, 
which tend to more resemble rounded cylinders or spheroids 
with roughly 4: 1 aspect ratio. Unlike for scalar wave problems, 
where exact solutions can be generated also for ellipsoidal 
targets, the vector field Maxwell equations fail to separate in 
ellipsoidal coordinates [3] and fully analytic solutions do not 
exist. 

As described in more detail below, the modeling approach 
applied here uses simplifications available for UXO-like target 
shapes, and also in different target electrodynamic regimes, to 
generate a combined prediction that quantitatively describes 
the full response. We specifically consider TDEM induction 
measurements. Here the transmitter loop current pulse gener- 
ates a magnetic field in the target region, and this changing 
applied field, especially as the pulse terminates, induces cur- 
rents in the target, generating a scattered magnetic field. The 
decaying scattered field, following pulse termination, induces 
the measured voltage in the receiver loop. 

In such a measurement there are three different regimes that 
one may identify in the voltage time traces: early, intermediate, 
and late time. At very early time, immediately following 
pulse termination, the currents are confined to the immediate 
surface of the target. The initial diffusion of these currents 
into the target interior leads to a power law decay (l/t^^"^ 
for nonferrous targets, for ferrous targets [4], [5]). At 

intermediate time, as the currents penetrate the deeper target 
interior, the power law crosses over to a multi-exponential 
decay, representing the simultaneous presence of a finite 
set of exponentially decaying modes [6]-[8]. Finally, at late 
time only the single, slowest decaying mode survives. We 
have developed a highly efficient combination of analytic 
and numerical models, based on rigorous solutions to the 
Maxwell equations, that covers these three regimes, and the 
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central purpose of this paper is to validate these models 
against laboratory data from both artificial spheroidal and real 
UXO targets, and to perform some inversion experiments that 
support their use for target discrimination. 

The outline of the paper is as follows. Details of the EM 
theory underlying the models, and their numerical implemen- 
tation, has been presented elsewhere [5], [8], but a basic 
overview is given in Sec. II. In Sec. HI the basic parameters 
of the NRL TEMTADS system are detailed. In Sec. ' model 
predictions are compared with TEMTADS data for spherical 
targets, for which an exact analytic theory also exists (Sec. 
IV-A), and for prolate (elongated) and oblate (discus-like) 
spheroidal targets (Sec. IV-B). In Sec. V, we describe results 
for the inverse problem, in which various target properties are 
treated as unknown, and seek to extract them from the data. 
In Sec. VI we describe results for certain real UXO targets 
(specifically, 60 mm and 81 mm mortar bodies). Finally, 
conclusions and directions for future work are presented in 
Sec. VII. 

II. Modeling Background 
A. Intermediate- to late-time modeling: mean field approach 

Our approach to the intermediate and late time regimes 
is based on a perturbation expansion about low frequency 
that takes advantage of the fact that analytic solutions for 
elUpsoidal targets do exist in the electrostatic limit (where 
the electric and magnetic fields are gradients of scalar fields). 
Based on this, we have developed a perturbation expansion 
about low frequency [6]-[8] that has an extremely efficient 
numerical implementation. The theory is dubbed the "mean 
field approach," since the expansion is highly nonlocal in 
space, with the currents and fields at any given point in the 
target being sensitive to their values throughout the target. 
Although formally valid only at low frequency, the theory is 
extended to higher frequencies by generating a large number 
of terms in the series (for a related numerical approach using 
an expansion in spheroidal wavefunctions, see also Refs. [9]- 
[11]). 

For time-domain measurements, low frequency corresponds 
to later time, in which initial rapid transients have died away. 
The solution to the Maxwell equations allows one to represent 
the electric field following pulse termination as a sum of 
exponentially decaying modes, 

oo 

E(x,t) = ^A„e(")(x)e-^"* (I) 

n=l 

where A„ are decay rates, e^"^ are mode shapes, and An are 
excitation coefficients. The first two are intrinsic properties 
of the target, analogous to vibration modes of a drumhead. 
Only the excitation amplitudes depend on the details of the 
measurement protocol. At early time a very large number of 
exponentials is present, and in fact the previously mentioned 
power laws arise from this large superposition (see Sec. II-B 
below). 

As time progresses, modes with larger values of A„ decay 
more quickly, and so at any given time t the signal will be 
dominated by some finite set of modes, namely those modes 



with A„ < 1/t. At very late time, t > 1/Xi, only the slowest 
decaying mode contributes, and the signal becomes a pure 
exponential decay. Thus, the earlier in time one wishes to 
model quantitatively, the greater the number of modes that 
are required. The ultimate limitation turns out to be the rate 
at which the excitation in pulse is terminated. If the pulse 
is turned off on a time scale tr (see Sec. lU-B), then only 
modes with < 1/ir have substantial amplitudes An, and 
a finite set of modes suffices for a full description of the 
target electrodynamics. For large targets, this may require 
many thousands, or even tens of thousands of modes, which 
is beyond current computational capability. 

However, for computational purposes it is only required 
that enough modes be computed that the resulting multi- 
exponential series overlaps the early-time regime. The early- 
time power law and mode descriptions may then be combined 
to fully describe the target dynamics over the full measured 
time range. For modes that decay slowly enough, hence 
contain low enough frequencies, the mean field approach can 
be used to compute them, and compute as well the excitation 
level of each. We will see that a few hundred modes is more 
than enough to attain the required overlap, and this basically 
serves to define the beginning of what we call the intermediate 
time regime. 

Using the mode orthogonality relation (which follows from 
the Maxwell equations), 

y"d3a;a(x)e(™)*(x) • e(")(x) = Smn, (2) 

where (t(x) is the conductivity, the excitation amplitude can 
be determined as 

An=4"^NT [ e(")*(x).dl, (3) 

in which the transmitter loop has been approximated by an 
ideal ID loop Ct with Nt windings, and 

4"^ = - /" dte^-'dtlT{t) (4) 

depends on the history transmitter loop current It (t) up until 
the beginning of the measurement window, taken here as t = 

0. To gain some intuition, a single perfect square wave pulse 
of amplitude and duration tp, one obtains 

4") = 7^(1-6-^"*-). (5) 

For a mode that decays rapidly on the scale tp, one has A„tp ^ 

1, and 7^"^ ~ 7^. For a more slowly decaying modes, 7^"-' will 
have a strong dependence on tp and n. In fact, for large targets 
one may actually encounter for small enough n the regime 
Xntp < 1 [e.g., tp = 25 ms and t„ = 1/A„ = 0(100 ms)] 
where 7^^^ will depend not only on tp, but on previous pulses. 

Finally, the measured voltage takes the form 

oo 

Vit) = E Ke-^"* (6) 

71 = 1 

in which, approximating the receiver as well by an ideal ID 
loop Cji with Nu windings, the voltage amplitudes are given 
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by the line integrals 

Vn=A,,NRf e(")(x).dl. (7) 

Equations (3)-(7) provide all the required ingredients for 
generating predicted data based on a target and measurement 
platform model. Our "mean field" numerical code divides 
naturally into two parts. 

The internal code solves the Maxwell equations to produce 
the intrinsic mode quantities A„ and e'"-' for a range of 
expected targets. With increasing A„, the modes have more 
complex spatial structure, and finite numerical precision means 
that only a finite set (a few hundred) of slowest decaying 
modes are actually produced [7], [8]. 

The external code uses the mode data, along with the 
measurement platform data, to compute current integrals ( ), 
the line integrals in (3) and (7), and then combines them to 
output the voltage amplitudes Vn and hence the time series (d). 
Note that the line integral computation requires full knowledge 
of the relative position and orientation of the target and 
platform. 

For high precision, the internal code can take anywhere from 
minutes to hours to produce mode data for a single target. 
However, given this data, the external code takes at most a 
few seconds to produce the full predictions. Precomputation 
and storage of a rapidly accessible database of target data is 
therefore essential. 

B. Complementary early time modeling 

For a rapidly terminated transmitter pulse, the external 
electric field, and induced voltage, display an early time power 
law divergence [4], [5] (saturating at very early time only 
on the scale of the off -ramp time t^. The boundary between 
the intermediate (multi-exponential) and late time (mono- 
exponential) regime occurs at the diffusion time scale 

TD = L^/D (8) 

where L is the characteristic target radius, and D = c^/ATTfia 
is the EM diffusion constant — this is the time scale required 
for the initial surface currents to diffuse into the center of the 
target. The early time regime corresponds to times t <C td 
(say, t < T£)/100), beginning deep into the multi-exponential 
regime where many (e.g., hundreds of) modes are excited. In 
this regime, for nonpermeable, or weakly permeable targets 
(yU ~ fib), one obtains the simple power law prediction 
prediction [4] 

Vit)=Ve/t^/^ t<^TD, (9) 

with all of the target and measurement parameters encom- 
passed by the single amplitude Ve, whose computation requires 
the solution of a certain Neumann problem for the Laplace 
equation in the space external to the target. 

For permeable targets, a new magnetic time scale 

Tmag = TnifJ-b/lJ-)'^ (10) 



emerges. For ferrous targets, /i//if, = O(IO^), and T-^-^ag/Tc ~ 
O(10~^) is tiny. The early time voltage then has a more 
complex magnetic surface mode structure, 

oo 

y{i) = Y.^nH{nr.Vt) (11) 

n=l 

where the k„ are surface mode eigenvalues, and the mode time 
trace profile 

H{s) = — ^ e*^crfc(s) 

^ I (12) 

where erfc(s) is the complementary error function, interpolates 
between a \/t^/'^ power law at early-early time, t ^ Tmag, and 
a power law at late-early time, T,nag <^ t <^ T£,. For 

large ferrous targets, this latter interval is very large, and may, 
in fact, accurately represent the signal over nearly the entire 
measurement interval (see Sec. IV). 

Figure I illustrates the important features of the early time 
modeling, including the complex evolution of the surface 
current depth profile [which extends H{s) to a function of both 
time and space [5]] that ultimately gives rise to the externally 
measured voltage (II). 

The surface modes are special surface current profiles (two 
such patterns are shown in Fig. 12 below) that, instead 
of decaying exponentially, evolve according to the universal 
function H{s). They and the k„ are solutions to an eigenvalue 
problem defined on the surface of the target [5]. They may be 
determined analytically only for spherical targets, where one 
finds 

Kl = W^mag, / = 1,2,3, (13) 

each (21 + 1) -degenerate, with Tmag = Anafila^ / fic'^, where 
a is the radius. The surface current patterns are controlled by 
the spherical harmonics of order /. The amplitudes V^'^ again 
require a solution to an external Laplace-Neumann problem. 

Unlike the bulk, exponential modes, under most conditions, 
only a very few surface modes are excited. The initial surface 
current pattern more-or-less follows the shape of the magnetic 
field generated by the transmitter coil. Unless the target is close 
to the coil, this field is fairly uniform, and the corresponding 
surface current density is fairly uniform as well, and can then 
be represented by the first few (two or three) modes. There is 
a very heavy numerical overhead in computing these modes 
and their excitation amplitudes, all in pursuit of predicting the 
rather limited information content of just a few coefficients. 
Given the success of extending the mean field predictions 
into the intermediate-early time regime, we have therefore 
found that it is much more efficient to extend the voltage 
curve by fitting the data at intermediate times to a one or 
two term series of the form (11), estimating k„ w 1/ ^Tmag 
for the first few modes. Although this precludes quantitative 
predictions at early-early time, it provides an enormously 
useful qualitative confirmation that the functional form H{s) 
accurately describes the data. 
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Fig. 1 . Illustration of the early time evolution of the surface density depth profile from the target surface for nonmagnetic (left) and magnetic (center) targets, 
beginning from a delta-function initial condition (perfect step function pulse termination). Distance R = r/L is scaled by the target size, time r = t/ro 
by the diffusion time, so that k here corresponds ftn\/Tc in (I I), and is essentially the permeability contrast {fi — fn,)//!},. The profiles are plotted for a 
sequence of 26 equally spaced scaled times 10^* < "J" < 0.05 (earlier times coiTesponding to narrower profiles). The nonmagnetic profile exhibits a pure 
Gaussian spreading into the target interior, while the magnetic profile is much more complex due to the surface magnetic boundary condition. Its maximum 
is pushed inwards from the boundary, and decays more rapidly with time. The right plot shows the time trace for the current density at the surface, R = 0, 
and is essentially the profile H{Ky/T), equation (12), which appears in the measured voltage (',]). For k = (solid blue line) the r-dependence follows an 
exact 1/\/t power law. For ft > (solid red line) the r-dependence crosses over from the identical 1/\/t form at early-early time to the l/r"^/^ power law 
(dashed red line) at late-early time [the asymptotic forms displayed in (12)]. 
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Fig. 2. Sketch of NRL TEMTADS array consisting of a 5 X 5 array of 25 
independent, concentric transmitter and receiver coils, numbered from to 
24 as shown. Due to rapid decay of signals with target depth, precise (cm 
level) geometry and placement of the coils (summarized in Table 1) can have 
significant effect on the overall measured voltage amplitude. 



Sensor center horizontal separation 40 cm 

Transmitter coil center height 4.3 cm 

Transmitter diameter 35 cm 

Number of transmitter coil windings Nq^ 35 

Receiver coil center height 0.4 cm 

Receiver diameter 25 cm 

Number of receiver coil windings Nji 16 

TABLE I 

nrl temtads array geometry. the transmitter coil windings 
are 7.8 cm tall with 0.4 cm thick endcaps on top and bottom. 
Height is measured from the bottom side of the lower endcap, 

and the transmitters are then modeled as an idealized id 
square loops at 0.4 + 3.9 = 4.3 cm height. the receiver coils are 
vertically compact and lie at the bases of the transmitter 
coils, hence are modeled as idealized id square loops at 0.4 
cm height. 



III. TEMTADS PLATFORM 

A. Platform geometry 

The 5 X 5 NRL TEMTADS sensor array is sketched in Fig. 
2, and its geometrical parameters are summarized in Table 1. 
The loops Gt and Cr are all modeled as perfect squares with 
35 cm and 25 cm edges, respectively. The origin is taken to 
be at the base of the lower endcap for sensor 12, the positive 
X-axis towards sensor 13, the positive y-axis towards sensor 
7, and the positive z-axis vertically upwards. The transmitter 
and receiver loop centers then all have x- and y-coordinates 
that are multiples of 40 cm. The transmitters are all at z = 4.3 
cm, and receivers are all at 2 = 0.4 cm. Target positions and 
orientations quoted in later sections are all defined relative to 
this frame of reference. 

The precise overall voltage amplitudes, required at least for 
initial verification of the instrument calibration, turn out to 
be surprisingly sensitive to small changes in these numbers. 
The scattered fields are approximately dipolar, and the voltage 
therefore decreases roughly as l/d^ with depth d. For example, 
therefore, a 1 cm error for a 30 cm deep target then leads to 
a 20% error in the voltage amplitude. A consistent systematic 
error of this magnitude, in fact, is what led us to discovering 
the existence of the endcaps, and the vertical offset between 
the transmitter and receiver loops. 

B. Transmitter waveform 

The TEMTADS bipolar pulse sequence is shown in Fig. 
.1. Each pulse is tp =- 25 ms long, followed by a 25 
ms measurement window. An adequate model of the pulse 
waveform is the form: 



m 



c(l 



(14) 



with exponential onset time constant ti = 0.33 ms, off-ramp 
time t,. ~ 10 /is, and current amplitude /max = 5.7 ± 0.3 
a. This form misses some detailed multi-exponential behavior 
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Fig. 3. TEMTADS transmitter cuiTent bipolar' pulse waveform. Top: multiple 
periods. Bottom: single 100 ms period. 

during the pulse onset that can be shown to have negligible 
effect on the excitation coefficients. The second half of the 
full bipolar pulse, beginning at t = 2tp, is the same as the one 
above, but inverted. The functional forms in (14) are simple 
enough that analytic forms for the current coefficients (-i) may 
be computed straightforwardly. 

IV. Data comparisons 
A. Spherical targets 

Having described the electromagnetic model, and the plat- 
form model required to implement it, we now turn to its 
validation with real data. We begin with spherical targets, for 
which exact analytic solutions exist in both the early time 
[4], [5] and multi-exponential regimes [2]. This allows one to 
validate the sensor model under conditions where the target 
model is fully specified. 

Figure shows results for a 15 cm diameter aluminum 
sphere, plotted on both linear and log time scales — the latter 
much more clearly verifies the asymptotic 1/Vt early time 
power law. The agreement is quite remarkable — note that the 
vertical scale is in millivolts, not an arbitrary scaled unit. The 
only real fitting parameter is the conductivity, and the chosen 
value (7 = 3 X 10^ S/m is well within the range expected 
for aluminum. As discussed in Sec. Ill, the overall pulse- 
to-pulse transmitter current amplitude is stable only at the 
10% level. This leads to an identical uncertainty in the overall 
voltage amplitude. In the figure, an overall factor of 1.03 has 
been applied to the data to obtain an optimal fit, well within 
this uncertainty. The slowest decaying mode for this target is 
Ti = 21.5 ms, so the measurement window here barely enters 
the late time regime t > Ti. The mean field prediction, based 
on an approximate calculation of the first 232 modes [7], [8], 
is seen to accurately describe the data well into the early time 
regime. 



Fig. 4. Data and theory for a 15 cm diameter aluminum sphere with center 
lying 16.5 cm below the center of sensor 12 (see Fig. 2), which is also the only 
active sensor. The solid red line is the data, the dashed red line the prediction 
from the exact analytic solution for the sphere, the dotted red hne is the mean 
field prediction (based on 232 modes), and the dashed black line is the early 
time l/Vt power law. The 1.03 overall multiplier listed in the legend has 
been apphed to the data to optimize the fit, and is well within the expected 
10% fluctuation in the TEMTADS current amplitude. The vertical dashed 
line marks the rough division between the early time and multi-exponential 
(< 100 modes) regimes, and it is seen that the mean field prediction is valid 
well into the early time regime. The slight deviation of the data from the 
analytic prediction at very eai'ly time, t < 0.1 ms, is likely an instrument 
saturation effect (seen much more clearly in Fig. 5, beginning roughly at the 
same voltage level). 



The mean field prediction has much more interesting struc- 
ture for ferrous targets. Due to the nature of the EM boundary 
conditions in the large permeability contrast limit, rather than 
computing only the slowest decaying modes, two distinct sets 
of slow (169 modes in this case, with time constants larger than 
3.01 ms) and fast (63 modes in this case, with time constants 
smaller than 0.74 ms) decaying modes are produced, with 
large gap between that would only be filled if one pushed the 
computation to higher order This is the source of the S-curve- 
like structure seen in the right panel of Fig. 5. The reduction 
in the number of slowly decaying modes reduces the accuracy 
of the theory near the early-intermediate time boundary (as 
compared to the nonmagnetic case shown in Fig. 4), but 
the presence of the more rapidly decaying modes at least 
provides an improved trend at very early time. The slowest 
decaying mode for this target has a time constant ri = 180 
ms, indicating a late time regime an order of magnitude beyond 
the measurement. 

The early time prediction, which follows both the exact 
solution and the data over a significant fraction of the time 
interval, deserves some comment. As described in Sec. 11-B, 
to obtain the solid black lines in Fig. ) the known eigenvalues 
(13) are used, but the amplitudes are determined (11) by 
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Fig. 5. Data and theory for a 15 cm diameter steel sphere with center lying 
16.5 cm below the center of sensor 12 (see Fig. 2), which is also the only 
active sensor The solid red line is the data (which shows clear instrument 
saturation above a few volts), the dashed red line the prediction from the 
exact analytic theory solution for the sphere, the dotted red line is the mean 
field prediction (based on 232 modes). The dashed black line is a two term 
fit to the early time form (II) using the known values (13), and the dotted 
black line is a single term fit using ki as a fit parameter. The vertical dashed 
line marks the rough division between the early time and multi-exponential 
(< f 00 modes) regimes, and is much later here than in Fig. ' because the EM 
time scale is proportional to the product crfi, which is an order of magnitude 
larger here. For reasons described in the text, the mean field prediction has 
a more complex structure for fen'ous targets, and penetrates only to the edge 
of the early time regime (it is the fact that it is accurate beyond about 20 
ms that is the real figure of merit here, as would be more evident if the data 
extended to later time). The 1.3 multiplier listed in the legend is that applied 
to the data to optimize the fit, and lies outside the expected 10% fluctuation in 
the cun'ent amplitude. The difference is likely the result of small positioning 
errors. Sensor saturation is apparent below about 0.5 ms. The late-early time 
l/t'"^^^ power law is evident in the data, but full convergence to the f/\/I 
early-early time power law is incomplete, and not expected until about 10 fis. 



fitting to the data. Only two terms are kept, 



Vit) = Vo 



( 



5) 



with the known value imag = 0.35 ms, and the amplitude 
Vo = 83 V, and mixing parameter a — 0.4 are fit. The one 
term series Vb = 60 V, a = provides an adequate, but lower 
quality fit. 

However, a better fit than both of these is provided by a 
single term series in which one allows the eigenvalue ki to 
be adjusted. The dotted black line in Fig. ) shows the result 
obtained using ki = l^s / \/*mag with l^s = 1.2, along with 
amplitude Vf = 80 V. This will be our fitting method of choice 
for non-spherical targets, where the eigenvalues k„ have not 
yet been computed. 

It is worth emphasizing the importance of the fact that 
analytic functional forms of the type (15) fit the data so well. 
The log-time plot demonstrates that the data span the full 
range over which the argument s in (12) interpolates between 
the two power laws [12]. The data therefore has significant 
structure through this time range, but this does not reflect any 
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Fig. 6. Consolidated plot of data and theory for a range of artificial aluminum 
prolate spheroidal targets, centered under Tx/Rxl 2 at various depths and 
orientations. The dimensions listed in the legend are diameters. Orientation 
angles indicate symmetry axis declination (towai'd the center of sensor 11), 
so that 0° con'esponds to vertical and 90° to horizontal. The multipliers are 
again the overall factors applied to the data to obtain optimal agreement with 
the prediction. The thick dashed lines are the mean field predictions, which 
show remarkable agreement well into the early time regime, where the onset 
of the power law is evident. 



deep structure of the target (beyond the fact that it is ferrous). 
Quite the contrary: as illustrated in Fig. : it represents the 
dynamics of a laterally very smooth surface current sheet as 
it begins to penetrate the first centimeter or so into target. 
The complexity arises strictly from the interplay between the 
electric and magnetic field boundary conditions at the surface. 
This serves to confirm that the early time regime provides 
limited target discrimination abihty (again, beyond the fact 
that it is ferrous). 

B. Prolate and oblate spheroidal targets 

Having verified instrument calibration and several other 
quantitative details under conditions where an exact solution 
exists, results for spheroidal targets are now presented. 

Figure 6 shows a consolidated plot of data and theory for 
various prolate (elongated) spheroidal aluminum targets at 
various depths and orientations. Spheroid aspect ratios Qz/axy 
vary between 2 and 5. 

The theoretical plots (thick dashed lines) are the mean field 
predictions based on the first 232 modes. It is again evident 
that the mean field predictions are valid well into the early time 
regime. The multi-exponential time series eventually saturates 
and falls below the data, but not before the l/Vt power law is 
quite well established. For smaller targets (e.g., the 4 x 4 x 20 
cm and 5 x 5 x 20 cm spheroids) the mean field prediction 
can cover nearly the entire measurement window. The multi- 
exponential time series eventually saturates and falls below 
the data, but not before the 1/v^ power law is quite weU 
established. Interpolating between the mean field prediction 
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Fig. 7. Consolidated plot of data and theory for a range of artificial 
aluminum oblate spheroidal targets, centered under Tx/Rxl2 at various depths 
and orientations. The thick dashed lines are the mean field predictions, which 
show remarkable agreement well into the early time regime, where the 1/Vt 
power law takes over (dotted lines). 



and this power law clearly enables one to accurately match 
the data over the full range (this will be demonstrated more 
quantitatively in the inversion experiments described below). 

It is apparent that most of the target discrimination infor- 
mation occurs at intermediate to late time. The traces are all 
more-or-less parallel at early time, and variations in the overall 
ampUtude from from variation in depth, size or geometry of 
the target are not distinguishable. On the other hand, at later 
time, the traces for smaller targets (e.g., again, the 4 x 4 x 20 
cm and 5 x 5 x 20 cm spheroids) drop off much more quickly 
than those of larger targets. 

There are also interesting dependencies on target orientation 
in this regime (green, red, magenta, and cyan curves in the 
upper part of the plot for the 10 x 10 x 20 cm spheroid 
[13]). For a vertical target, the excited modes are dominated 
by currents that circulate around the symmetry axis, while for 
a horizontal target the currents tend to circulate along it. The 
horizontal target mode has a slower decay rate (time constant 
Til = 13.7 ms vs. Ty = 12.0 ms), and couples differently to 
the transmitted field, and this is visible in the later-time traces. 

Identical conclusions are evident from the data on oblate 
(discus-like) spheroidal aluminum targets (aspect ratios a ~ 
0.2,0.4) shown in Fig. 7. Here we have overlayed segments 
of 1/Vi power law on each curve, explicitly demonstrating 
successful interpolation (with, perhaps, 5-10% errors in the 
overlap regime). 

The dependence on orientation is much stronger for oblate 
spheroids (green, red, magenta, and cyan curves for the 
20 x 20 x 8 cm spheroid [13]). Because it is being "squeezed" 
vertically, the horizontal target (discus on edge) mode now 
has significantly faster decay rate than vertical target (discus 
lying flat) mode (time constant Th = 13 ms vs. Ty = 24 ms). 
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Fig. 8. Consolidated plot of data and theory for a range of steel prolate 
spheroidal targets, centered under Tx/Rxl2 at various depths and orientations. 
The thick dashed lines are the early time predictions, which show remarkable 
agreement over nearly the entire measurement window. The latter take the 
form (11) with a single term, in which the amplitude Vj^ and eigenvalue 
'^1 = 'eff/Vmag are adjusted to optimize the fit. Here Tmag is defined 
by (10) and (x), with the choice L = mm{axy,az} = a^y. The mean 
field predictions are shown by the dotted lines. If extended over the full time 
interval, they also would display the S-curve behavior seen in Fig. 5. For these 
smaller targets, their region of validity begins only at later times where the 
signal levels are falling into the noise floor. The early time regime therefore 
encompasses almost the full range of useful data. 



Tx1 2, Rxl 2: St oblate spheroids under Tx/Rxl 2 at various depths-to-center 
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Fig. 9. Consolidated plot of data and theory for a range of steel oblate 
spheroidal targets, centered under Tx/Rxl2 at various depths and orientations. 
The thick dashed lines are the early time predictions, which show remarkable 
agreement over nearly the entire measurement window. The latter take the 
form ( i 1 ) with a single term, in which the amplitude Vj^ and eigenvalue 
Kl = leff/ y/Tuiag are adjusted to optimize the fit. Here Tmag is defined by 
(10) and (8), with the choice L = mm{axy , az} = a^. The mean field 
predictions are shown by the dotted lines. For these smaller targets, their 
region of validity again begins only as the signal levels are falling into the 
noise floor The early time regime (t < to-m) encompasses nearly the full 
range of useful data. 
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Because the the latter mode is not excited at all when the target 
is horizontal, the 90° (cyan) curve in Fig. 7 decays much more 
quickly at late time than the other curves. 

In both Figs. 6 and 7 the multiphers used to scale the data 
for optimal fit appear to have a small 10%) systematic 
bias that cannot be explained by random variation in the 
transmitter loop current. A combination of small conductivity 
and positioning errors is the likely culprit. 

Figures ^> and v show data and theory for steel prolate and 
oblate spheroidal targets. As for the steel sphere (Fig. 5), the 
early time regime dominates, and the mean field results (dotted 
curves; with S-curve behavior excised in this case so as not to 
confuse the plots) are valid only over a small part of the time 
interval where the data is already becoming quite noisy. In 
most cases, however, the fact that the data is dropping below 
the early time curve is evident, pointing to the necessity of 
a multi-exponential description. As before, these predictions 
actually push quite deeply into the early time regime, but 
the measurement window, and instrument dynamic range, are 
such as to strongly limit the information content of the multi- 
exponential part of the signal. 

V. Inverse problem for spheroidal targets 

We have so far described applications of the early time and 
exponential models to the forward problem, in which detailed 
model predictions for a known target are compared to data. 
Having demonstrated the quantitative success of the models 
for this problem, we now turn to the inverse problem, in 
which some set of target characteristics is treated as unknown, 
and one attempts to determine them by searching for the 
target model whose predictions best match the data. We have 
previously presented solutions to the inverse problem using 
noise corrupted simulated data [6]. Here we will base the 
inversions on the TEMTADS data. 

Of particular interest are ambiguities in the data, i.e., target 
properties that are poorly constrained by a particular data 
set, either due to poor data quality (e.g., noise), or due to 
fundamental trade-offs between certain parameters that exist 
even for essentially perfect data. We will see, for example, that 
it is very difficult to simultaneously determine target depth, 
size, and conductivity. We will also see that the enormous 
differences seen between magnetic and nonmagnetic target 
data in Sec. IV give rise to similar differences for the inverse 
problem. 

A. Objective function 

Inverse problems are generally formulated as the solution 
to the following optimization problem. Let m e M be a 
vector of forward model parameters spanning a model space 
M, in our case the space of target and sensor platform 
parameters that define the forward problem. Let d be the 
measured data vector, in our case the set of voltage time series 
corresponding to a given set transmitter-receiver combinations 
and platform positions. Let dpi.cd(m) be the forward model 
prediction for the data given a model m. For a perfect model 
and perfect data, there should be a unique model mo for 
which dpicd(mo) = d. In the presence of noise, and/or an 



incomplete model (e.g., real UXO are not ideal spheroids), 
such an exact match is not achievable, and we instead seek an 
optimal solution rather than a perfect solution. 

The sense in which a solution is optimal is defined by an 
objective function i^(m|d), which, e.g., vanishes if dpiod = d, 
but more generally attains a minimum value mo in the space 
of available models: 



mo 



arg min F(m|d). 



(16) 



Finding this minimum entails a search over the space M, and 
there may be multiple local minima that interfere with the 
discovery of the absolute minimum. It is generally the case, 
as well, that the problem is under-determined, and an entire 
subspace of different m, for example, may very nearly achieve 
the same minimum. To some degree, such ambiguities can be 
controlled using a priori constraints (e.g., on a target's per- 
mitted range of depths or conductivities), effectively biasing 
F(m|d) toward a preferred set of solutions. The danger, of 
course, is that this blinds one to targets disagreeing with this 
bias, and mo may end up being vastly different from the truth. 
We will see examples of this below. 

In what follows, we will use the following form of objective 
function: 



i^(m|d) 



Nt 

E 

k=l 



Wk 



Vprod(ifc;m) 



1 



Fp,ior(m), (17) 



in which tk, k = 1,2, Nt are the measurement platform 
time gates, Vmeas, ^prod are the measured and predicted 
voltages, and Fpiior contains any prior constraints. The weights 
Wk may be used to fine tune the weighting given to data 
in the different time regimes, for example suppressing noise- 
dominated time gates. The ratio l^neas/V^icd is used (adopting, 
basically, a logarithmic rather than a linear voltage scale), 
in place, e.g., of a more conventional mean square error 



in order to give more democratic weight 



to early and later-time regimes. Specifically, signals may be 
orders of magnitude weaker at later time, but it will be 
seen that the multi-exponential decay still contains key target 
identification information not present at early time. In addition, 
l^prcd is placed in the denominator because noise effects may 
easily induce sign changes in Kncas at later time, which would 
produce singularities in Vprcd/Kncas- 

B. Inversions for aluminum targets 

We begin with a set of numerical inversion experiments 
using the 10 x 10 x 20 cm aluminum spheroid data (up- 
per curves in Fig. 6). We consider some interesting issues 
involving tradeoffs between target depth, conductivity, and 
orientation, which are most clearly elucidated by treating the 
target geometry (in this case, axy ~ 5 cm radius, aspect ratio 
a = az/a^y = 2) as known. We have performed inversions in 
which the target geometry also varies (see also Ref. [6]) and 
found similar effects. We will also use the (highest quality) 
data only from the pair Txl2-Rxl2. In Sec. VI, some results 
will be shown using multiple Tx/Rx pairs. The inversion code 
uses the standard simplex method (which 'walks' its way 
toward a local minimum, in sequentially decreasing steps). 
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Fig. 10. Inversion experiment using TEMTADS data for a 10 X 10 X 20 
cm prolate spheroid, using the objective function (17). Prior information is 
incorporated in that only the parameters m = {rr, ri, te-m} ai'e permitted 
to vary. Target geometry and orientation (symmetry axis vertical) are fixed, 
and the transmitter current amplitude is fixed at 5.5 a. The upper plot shows 
the data (thick black line — the same as the green curve in Fig. d) and the 
remarkably accurate optimal fit (solid red line to the right of te-m; dashed 
red line to its left). The lower plot shows the convergence of F, cr, d with 
iteration number. Note the two-stage convergence, in which d equilibrates 
first, followed by the 'weaker' parameter cr. 



that does not require any gradients of the objective function. 
More sophisticated search methods, that may operate more 
efficiently, will be left for future work. 

The transition between the early-time and multi-exponential 
regimes for non-magnetic targets is treated as follows. One of 
the inversion parameters is taken to be a crossover time tc-m> 
and for a given value, the predicted voltage takes the form 

-1/2 ^ ^ ^ 



t < u. 



VMF{tc-m){t/tc-m) 

in which VuFit) is the direct mean field prediction based on 
the remaining parameters in m [14]. 

In a number of numerical experiments, it was found that if 
the target conductivity a, the depth d, and the tilt angle 6 are all 
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Fig. 11. Second inversion experiment using TEMTADS data for a 10 X 10 X 
20 cm prolate spheroid, using the objective function (17). The pai'ameters 
m = {d,9,ts-m} are now permitted to vary. Target geometry is fixed, the 
transmitter current amplitude is fixed at 5.5 a, and the conductivity is fixed 
at the value found in the first experiment (Fig. in). The upper plot shows the 
data (thick black line — the same as the magenta curve in Fig. h) and optimal 
fit (solid red line to the right of te-m; dashed red line to its left). The lower 
plot shows the convergence of F,9,d with iteration number. The value of 
^min = 63.4° is very close to the 60° ground truth. One sees even more 
clearly in this figure the initial equilibration of the depth, followed later by the 
convergence of the tilt angle, which has a much weaker effect on the voltage 
time series, mostly at later time. 



permitted to vary, the inversion is very unstable. Specifically, 
varying the tilt changes the signal decay rate at later time 
(reflecting differences in the decay rates of the excited modes). 
However, so does varying the conductivity, and the two effects 
can very nearly be made to cancel, so long as the depth 
is adjusted slightly to maintain the observed overall signal 
amplitude. The result is a solution mo with unphysical values 
of all three parameters. In realistic scenarios, one would then 
have to include in Fpiior terms that constrain, for example, the 
conductivity to in a range of acceptable values for aluminum. 
Since, for the data dealt with in this section, we have 
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B 



Fig. 12. Early time surface current flows depending on relative orientation 
of the magnetic field for a 4:1 prolate spheroid. Plotted is the stream function 
ip associated with two of the early time surface modes whose level curves 
are the stream Unes of the surface current (circulating clockwise around blue 
patches, counterclockwise ai'ound red patches). The characteristic radius reft 
associated with a mode is essentially the (half) distance between extrema of 
ip (red and blue regions). The diffusion of the surface current sheet into the 
tai'get, at any given point on the surface, is shown in Fig. I for different values 
of fc = (/i/47r(T)^''^creff . 

complete ground truth, we proceed in a slightly different 
fashion. We begin by determining an optimal value of the 
conductivity from one of the data sets using the known value 
of the tilt angle (as well as the target horizontal position). The 
depth and crossover time are permitted to vary as well. The 
results are shown in Fig. 10, where a value a = 2.84 x 10^ S/m 
is found. Using a different value of a is found to significantly 
change the inversion result for the tilt angle. 

Note that convergence of the inversion takes place in two 
stages. The depth equilibrates first, due to the strong 
dipolar dependence of the overall voltage amplitude. Once a 
reasonable value for d emerges, the more subtle conductivity- 
induced changes in the voltage curve shape (mainly at later 
time) can be productively optimized. 

Note also that the inverted depth dmin = 21.9 cm differs 
slightly from the 21 cm ground truth value. This may partially 
be due to experimental error, but is likely also due to the ~10% 
uncertainty in the transmitter current amplitude /max (see Sec. 
111-B) which has been fixed at the value 5.5 a. Inversion 
experiments have been performed in which /,nax is also treated 
as a free parameter, but this is also found to be highly unstable. 
Specifically, the voltage signal changes very little if /max 
is varied, while at the same time adjusting the target depth 
d — ^both mainly affect the overall amplitude of the voltage 
time series, not its shape [15]. For the same reason, here 
and in experiments described below, this ambiguity has very 
little effect on the inverted values of most other parameters, 
specifically those that indeed impact the shape of the time 
series. 

With an accurate conductivity value now in hand, in our 
second inversion experiment, we fix the former and invert for 
the tilt angle [16]. As seen in Fig. 1 1, this indeed produces an 
inverted f?min very close to the 60° truth. The two-step nature 
of the convergence is seen here as well. Similar experiments 
(not shown) using the 30° and 90° tilt data (red and cyan 
curves, respectively, in Fig. 6) produce equivalent results. The 
inverted depth dmin = 25.2 cm differs only slightly from the 
24.4 cm ground truth. 



We have performed inversions as well using the 20 x 20 x 8 
cm oblate spheroid data shown in Fig. '') which yield very 
similar results (not shown). One again fits the conductivity 
(and depth) using the 0° tilt data (green curve in Fig. 7), and 
then uses this value to invert for tilt angle (and depth) using 
the other data sets. One again finds delicate features where a 
slightly incorrect conductivity value leads to a highly inaccu- 
rate tilt angle. The problem actually appears to be somewhat 
worse in this case because of the early time crossover time 
te-m occurs relatively earlier for oblate geometries, increasing 
the model misfit for the given number of modes (232). 

C. Inversions for steel targets 

We next present inversions based on the data shown in Fig. 
9 for steel oblate spheroidal targets. These turn out to serve 
as an illustration of another set of inversion pitfalls one may 
encounter, in this case depending on how one balances prior 
knowledge in the presence of noisy data in the later time 
regime. 

Since the early time behavior of ferrous targets is much 
more complicated than that for nonmagnetic targets, in addi- 
tion to the crossover time tc-m, one must provide an appropri- 
ate parametrization of the magnetic surface mode series (11) 
and (12) [5]. As described in Sec. II-B, replacing the series 
by a single term with effective parameters is found to provide 
the best fit. For reasons that will become evident below, we 
parameterize it in the form 

Ve{t)=Vo'H{^t/7,), (19) 

where 



with Tcff an effective early time target radius. This parameter, 
derived as previously from the early time eigenvalue k„, 
reflects the primary length scale of the surface currents for 
a given mode, and decreases as n increases and the geometric 
complexity of the mode increases. This length scale may be 
thought of as the diameter of the target along the direction of 
the magnetic field that primarily serves to excite the mode. For 
example, currents circulating around the symmetry axis Qz are 
generated by a vertical magnetic field (left panel of Fig. 1 2), 
and Toff is then found to be comparable to a^. For magnetic 
field orthogonal to Uz the currents circulate up one side of the 
symmetry axis and down the other (right panel of Fig. 12), 
and Toff is found to be comparable to the radius axy- 

As a consequence the length Uz will be reflected in the early 
time data when Uxy is what would observed in a corresponding 
visual inspection of the target from the point of view of a 
transmitter-receiver above. Correspondingly, when the target 
is laid with symmetry axis horizontal, a^y will be exhibited 
when a visual inspection would clearly see both a^y and a^. 
These observations will be confirmed by the data below. 

Analogous to the inversion shown in Fig. 10, in Fig. 13 
we show inversions for the parameters d, a, tc-m, ''cff with tilt 
angle fixed at the 0° ground truth value, and fi = 100. As 
can be seen, the crossover time t^-m — 10 ms lies in a regime 
where noise effects are becoming significant. In particular, in 



11 



WEICHMAN: VALIDATION OF ADVANCED EM MODELS FOR UXO DISCRIMINATION 




Inversion convergence: 20x 20 x 4 cm St spheroid; d = 52 cm:9 = 0° 



■S 10' 



10" 



o 



100 150 200 250 300 

iteration number 



1 son 



S 45 




d = 43.7 cm 



50 



100 150 200 250 300 

iteration number 



o , = 0.704 X 1 0° SIrf 



c 50 100 150 200 250 300 

o iteration number 

(/} 

I' 
I" 

§ 50 100 150 200 250 300 

iteration number 



10 



Inversion convergence: 20x 20 x 4 cm St spiieroid; d = 52 cm; 9 = 0° 







F . = 1.9865 




mm 



52 
£■ 51 

£ 50k« _"%* 
S 49- 



50 



I 5 
S 4.5 

E 4- 

I 3-5 

S 3 - 







100 150 
Iteration number 



d . = 51 .1 cm 



100 150 
Iteration number 



r . = 3.35cm 



100 150 
Iteration number 



Fig. 13. Two inversion experiments using TEMTADS data for a 20 X 20 X 4 cm steel oblate spheroid, using the objective function (17). Rather than 
producing a simple power law, the early time theory now produces the nontrivial 'error function' series (II) and ( 1 2), modeled here as a single term with 
optimized value of the effective radius r^ff (solid red line in the upper plots). The target orientation is held fixed (symmetry axis vertical). Left hand plots: 
the depth, conductivity, and effective radius are permitted to vary. The resulting over-tit of the later-time noise produces unphysical values of all three. Right 
hand plots: the conductivity is now fixed at a physically reasonable value cr = 3 X 10® S/m. The depth, conductivity, and effective radius now converge to 
physically reasonable values as well. 



this case the noise happens to induce a sharper downturn in 
the signal in the multi-exponential regime than is case for the 
true signal. In attempt to fit this, the inversion on the left 
unphysically suppresses the conductivity (o-min = 7.04 x 10^ 
S/m), and compensates with an unphysically large effective 
radius (r,nin = 7.37 cm, much larger than the 2 cm half- 
thickness of the discus). This keeps the product ar^g in (Hi) 
essentially fixed, thereby maintaining a good fit in the early 
time regime. Note that the depth d^wi = 43.7 cm also differs 
substantially from the 52 cm ground truth. 



In order to avoid this overfit of the noise, the right hand 
side of Fig. 13 shows the superior result obtained by fixing 
the conductivity at the physically reasonable value cr = 3x 10^ 
S/m, with the results rmin = 3.35 cm and dmin = 51.1 cm. 

The main lesson here is that significantly different infor- 
mation (specifically, different combinations of conductivity 
with other parameters) is contained in the early time and 
later time regimes for ferrous targets, and there are large 
parameter ambiguities in the absence of good data in both. For 
ferrous targets, the more rapid 1 / 1^^^ dominant decay in early 
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Fig. 14. TEMTADS data (solid lines) and theoretical fits (dashed lines) taken for a 60 mm (left) and 81 mm (right) mortar body targets lying flat (red and 
magenta curves) and nose-up vertical (blue and cyan curves). The early time fits, based on ( l^>) and (20) with fj, = 100 and cr = 5 X 10^ S/m, produce the 
very reasonable effective diameters d^ft = 2rcS (see the illustrations in Fig. ! 2) as indicated in the legends. The length listed at the beginning of each legend 
entry is the target depth (ground truth for the data and early time curves; fitted value for the mean field curves). The mean field predictions (solid magenta 
and cyan curves), using these same parameters, properly extend the early time curves but completely fail to fit the data. Lowering the conductivity by more 
than an order of magnitude accurately mimics the increased rate of decay (dashed magenta and cyan curves), but the overall voltage amplitude is incorrect 
(as indicated by the applied multipliers in the legend). As described in the text, and in Fig. 15, the correct explanation for the discrepancy is the finite (0.5-1 
cm) mortar shell thickness. 



time leads to much smaller signals, hence degraded data, in 
the multi-exponential regime. Parameters relying on the latter 
will therefore be more poorly determined, and one may be 
forced to apply a larger set of prior knowledge constraints than 
for nonmagnetic targets. Since real UXO targets are usually 
ferrous, this lesson has important practical implications. 

VI. Some results for Real UXO targets 

We turn finally to some initial inversion results for real 
UXO targets, namely 60 mm and 81 mm mortar bodies, 
selected for their near-spheroidal shape. The former has a 6 
cm diameter, and is approximately 13 cm long; the aspect 
ratio is therefore taken as a = 2. The latter has an 8.1 cm 
diameter and is approximately 25 cm long; the aspect ratio is 
therefore taken as a = 3. Real UXO are not ideal spheroids, 
but we will see that key target discrimination information can 
be obtained by comparing their electromagnetic responses to 
those of spheroids with similar geometry. 

Data and theoretical fits for both UXO are shown in Fig. 14. 
The most important observation is that, with the assumptions 
(7 = 5 X 10^ S/m and fj, = 100, the early time curves 
for the two orientations (vertical and flat) produce excellent 
fits, and estimates for the effective radius that accord very 
nicely with the discussion in Sec. V-C (and Fig. 12): for the 
vertically oriented mortars, = 10.5, 14.7 cm are indeed 
comparable to the mortar half-lengths, while for 'flat' mortars, 
Tcff = 3.7, 3.9 cm are comparable to their radii. Note that 
these inferences are independent of the target depth (which, 
to leading order, affects only the amplitude Vq). 

The next observation is that, contrary to the early time 
fit, the mean field predictions (dashed blue and red lines 



in Fig. ) fail to fit the data. Thus, although they provide 
a logical extension of the early time curves [17], the data 
follow a much more steeply decaying path. As indicated by 
the red and magenta dashed lines, the decay can be fit using 
much smaller (by factors of 20-50) conductivities. However, 
these unphysically small effective values are also unphysically 
orientation dependent, and give similarly inconsistent values 
for the overall voltage scales (hence the x4-x7 fudge factor 
multipliers listed in the figure legends). 

The origin of these modeling discrepancies is that the 
finite mortar shell thickness (in the 0.5-1 cm range) has 
not been accounted for. To illustrate this, in Fig. 15 exact 
analytic results for a hollow sphere are shown. The results 
there confirm the rapid increase in the mode decay rates 
with decreasing shell thickness (left panel), and the resulting 
early breakaway of the voltage curves from the early time 
result (right panel). The breakaway time scales with the shell 
thickness as tx ~ ^shcii/-^' where D ~ c^/A-K^a is the 
electromagnetic diffusion constant (see Sec. II-B). With the 
given values, this relation may be put in the form; 

1 ms VI. 3 mm/ 
In Fig. 14, ones observes tx ~ 2-3 ms, which indeed places 
Wsheii in the right range (and is consistent as well with the 
ah — 0.8, 0.85 curves — hence lUshcii = 0.3-0.5 cm — in Fig. 
15). 

To summarize, the early time fits provide direct estimates of 
the target geometry, while the breakaway time tx away from 
the early time extrapolation (as well as the solid spheroid mean 
field prediction), provide direct estimates of the UXO shell 
thickness. These are key parameters in target identification. 
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Tx/Rx12: 6 cm diameter hollow steel spheres (d = 30 cm) 



E 





c; = 5x10®S/m 




^,, = 100 
rei 


solid sphere, = 




Early time (r^^^ = 2.5 cm) 






a^ = 0.6 




a, = 0.7 




a^ = 0.8 








= 0.85 




a^ = 0.9 




10-^ 10° 


10^ 



Time (ms) 

Fig. 15. Exact analytic TDEM voltage predictions for a sequence of hollow 
steel spheres, including the solid sphere (d^ = 0, solid black line). Since the 
surface geometry of all targets is identical, the same early time curve (dashed 
line) fits all data sets. However, the multi-exponential regime begins earlier 
for thinner-shelled tai'gets. 

Quantifying the latter would require an enhancement of the 
mean field code to deal explicitly with hollow spheroids. This 
can be done with modest effort, and will be considered for 
future work. 

A. Information content of multiple Tx-Rx combinations 

As described in the previous subsection, the early time 
prediction provides an effective orientation-dependent target 
radius. For a given buried target, the orientation is fixed, and 
one seeks other ways of extracting this information. Here one 
may take advantage of the array degrees of freedom available 
with the TEMTADS platform; each Tx-Rx pair effectively 
provides a different 'look' at the target. In Fig. 16 we show 
such data for 20 x 20 x 8 cm oblate steel spheroid (left panel), 
and the 81 mm mortar (right panel). In each case the targets 
are centered (50 cm and 48 cm, respectively) below Rx/Txl2, 
with their symmetry axes vertical. See Fig. 2 for the labeling. 

Note that if the transmitter and receiver coils were identical 
then, by complementarity, the Tx(m)/Rx(n) response would 
coincide with that of Tx(n)/Rx(m). The fact that they are 
somewhat different (see Table ) explains, for example, the 
small differences between the sets of green and blue curves 
(e.g., Tx7/Rxl2 vs. Txl2/Rx7 or Txl7/Rxl2 vs. Txl2/Rxl7). 
On the other hand, if two sets of Rx/Tx pairs are symmetrically 
placed relative to the (axially symmetric) target, then the 
responses should be identical (e.g., Tx7/Rxl2 vs. Txll/Rxl2 
or Txll/Rxl3 vs. Tx7/Rxl7), and this is indeed observed. 

To leading order, the actual time series is a superposition of 
the horizontal and vertical early time magnetic surface modes 
(Fig. 12), and one expects that a fit to a single value of 
r^g will then find values that interpolate between the two. 
This is indeed, for the most part, seen to be the case. For 
the combination Rxl2/Txl2 (red lines), the effective radius is 
indeed comparable to the 'vertical half-height' of the target 



(4 cm, vs. Toff = 4.3 cm for the oblate spheroid; 12 cm vs. 
Toff = 14.7 cm for the 81 mm mortar). 

So long as either Rxl2 or Txl2 is used (e.g., green and blue 
curves in Fig. 16), r^s remains comparable to (or even larger 
than) its Rxl2/Txl2 value. For Txl2, this is because the same 
dominant magnetic surface mode is being excited, irrespective 
of which receiver is used to observe it. For Rxl2 the same 
argument proceeds through complementarity: an off-center 
transmitter will excite more than one surface mode, but the 
symmetrically placed receiver will see only the symmetrically 
circulating mode. 

However, when both transmitter and receiver are off-center 
(e.g., cyan and magenta curves in the figure), the effective ra- 
dius is seen to become comparable to the physical (horizontal) 
radius of the target (10 cm for the oblate spheroid; 4.05 cm 
for the mortar). Thus, despite the fact that the overall signal 
levels drop precipitously as the Tx-Rx separation increases, 
to the point where the multi-exponential regime becomes 
unobservable [18], robust target geometry inferences can still 
be made using this data through the early-time prediction fits. 
It should also be noted that these fits, being restricted to the 
early time regime, entail adjustment of the single parameter 
Toff (or k), and so do not require a sophisticated inversion 
scheme. 

There is one feature of the 81 mm mortar data that deserves 
further comment. The Txll/Rxl3 and Tx7/Rxl7 (cyan curves 
in the left panel of Fig. ) display an unusual cusp feature 
at about 0.3 ms. Since the magnitude of the voltage is being 
plotted, this actually corresponds to a node in the response, 
i.e., a sign reversal of V{t). The easiest way to understand this 
effect is to note that a target close to Txl2 will generate a net 
downward-pointing magnetic field through Rxl3. However, 
as the target depth increases, the field becomes net upward- 
pointing, reversing the sign of the flux, and hence of the 
induced voltage. In this picture, a node in the response occurs 
as a function of target depth, but a similar argument would 
also produce a node as a function of time due to exchange 
of dominance of two early-time mode contributions. Vertical 
variation of the applied magnetic field will always produce 
higher order modes (with more complex spatial structure than 
those shown in Fig. I '') with faster decay rates, and near the 
critical depth the multiple contributions sum to produce a zero 
crossing. 

VII. Summary and conclusions 

The results presented in this paper demonstrate the unprece- 
dented accuracy available from our first principles, physics 
based models covering the entire measurement window, from 
the early time multi-power law regime, all the way through 
the multi-exponential regime to the late time mono-exponential 
regime. Prior to the mean field code's current iteration [7], [8], 
the number of accurately computed modes used to describe the 
multi-exponential regime was limited to perhaps a few dozen 
[6] . As seen in the results presented, by generating the required 
overlap of the early time and multi-exponential regimes, this 
improvement is critical to the success of the validation and 
inversion tests. 
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20x20x8 cm St olbate spheroid: 'flat' 50 cm under Tx/Rx12 



81 mm mortar: nose up 48 cm under Tx/Rx12 
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Fig. 16. Data (thin solid and dashed lines) and early-time model fits (thick solid fines) for two fixed steel targets using different Tx/Rx combinations (with 
labeling shown in Fig. ''). Left: 20 X 20 X 8 cm oblate spheroid. Right: 81 mm mortar body. For these symmetrically placed targets, the solid and dashed 
curve pairs of each color (e.g., Tx7/Rxl2 and Txll/Rxl2; Txll/Rxl3 and Tx7/Rxl7) should be degenerate, and this is indeed observed. In both cases the 
early-time effective target radius (listed in the legend) is seen to vary with 'look angle' in a way consistent with the target and measurement geometry. Note 
that the breakaway time (21) continues to be observable in a number of the 81 mm voltage traces (and continues to be notably absent in the solid spheroid 
traces). The cusp-like feature in the 81 mm Txll/Rxl3 and Tx7/Rxl7 traces (sofid and dashed cyan curves) actually represents a sign change in the voltage, 
as discussed in the text. 



It should be emphasized that the increase in predictive 
power continues to operate with extremely high numerical 
efficiency. The creation of the mode data for a given target 
cannot be performed in real time, but once this data is made 
available in a database that spans the expected target geome- 
tries, its acquisition and use for measurement predictions can 
be performed in real time — operating at essentially the same 
speed as predictions using the exact solution for the sphere. 

As seen in the figures, the dominant regimes visible in the 
data depend very strongly on the target size and physical 
properties. Increasing target size and magnetic permeability 
expands the early time regime to later physical time. Smaller 
aluminum targets (e.g., blue lines in Fig. 6) are completely 
described by the mean field approach over the full time 
range, while even the smaller steel targets barely enter multi- 
exponential regime (see Figs. 8 and 9) before the signal fades 
into the noise floor [19]. 

Through investigation of inversions it has been shown 
that different time regimes and different Tx/Rx pairs contain 
complementary target discrimination information. The most 
exciting development is the observation that the early time 
regime, which dominates the steel target response, contains 
direct information about the target geometry (via the different 
effective target diameters seen from different 'look angles') 
and the hollow target shell thickness (through the earlier 
breakaway time to multi-exponential decay for thinner shelled 
targets). These are key features that will support target iden- 
tification and clutter rejection. 

All of these features, whose quantitative interpretation is 
enabled by the present models, can be applied to the pur- 
suit of robust target discrimination and identification under 
more challenging conditions. The code efficiency becomes 
especially critical for this purpose, as searches through the 
database for the target whose response best matches the data 



requires hundreds, or perhaps even thousands, of iterations. 
In the inversion results presented, these searches presently 
take several minutes. Further focused improvements in code 
efficiency could probably reduce this to under one minute. 
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